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Abstract 

A rare human rotavirus, G3P[9] strain RVA/Human-tc/KOR/CAUl 2-2-5 1/2013/G3P[9], was isolated from tlie stool of a 9-year- 
old female hospitalized with acute watery diarrhea in August 2012 in South Korea using a cell culture system, and its 
genome was analyzed. The complete genomic constellation of the CAU 12-2-51 strain revealed a novel genotype 
constellation for human rotavirus, G3-P[9]-I2-R2-C2-I\/12-A3-N2-T3-E3-H3. Phylogenetic analysis revealed that the CAU12-2-51 
strain originated from feline- and bovine-like reassortment strains. The genes encoding VP4, VP7, NSP1, NSP3, NSP4, and 
NSP5 were related to human/feline-like and feline rotavirus strains, whereas the remaining five genes encoding VP1, VP2, 
VP3, VP6, and NSP2 were related to the human/bovine-like and bovine rotavirus strains. This novel strain was identified for 
the first time, providing evidence of feline/bovine-to-human transmission of rotavirus. The data presented herein provide 
information regarding rotavirus diversity and evolution. 

Citation: Jeong S, Than VT, Lim I, Kim W {201 4) Whole-Genome Analysis of a Rare Human Korean G3P[9] Rotavirus Strain Suggests a Complex Evolutionary Origin 
Potentially Involving Reassortment Events between Feline and Bovine Rotaviruses. PLoS ONE 9(5): e97127. doi:10.1371/journal.pone.0097127 

Editor: Karol Sestak, Tulane University, United States of America 

Received February 21, 2014; Accepted April 14, 2014; Published May 12, 2014 

Copyright: © 2014 Jeong et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) 
{2013R1A1 A2A10012148) and the Brain Korea 21 PLUS Program funded by the Ministry of Education. The funders had no role in study design, data collection 
and analysis, decision to publish, or preparation of the manuscript. 

Competing interests: The authors have declared that no competing interests exist. 

* E-mail: kimwy@cau.ac.kr 



Introduction 

Rotaviruses are the major cause of acute gastroenteritis in 
young chUdren worldwide. Each year, these viruses are associated 
with deaths of approximately 453,000 children younger than 5 
years of age, most of whom reside in low-income countries in 
Africa and Asia, where access to potable water, sanitation, and 
medical care is often limited [1]. Rotavirus, a member of the family 
Reoviridae, is classified into eight species/groups (A-H) [2]. Its 
genome is composed of 1 1 gene segments of double-stranded RNA 
that encode six structural viral proteins (VP1-VP4, VP6, and VP7) 
and six nonstructural proteins (NSP1-NSP6) [3]. 

The outer capsid of rotaviruses contains two proteins, the 
glycoprotein VP7 and the protease-sensitive VP4, which are used 
to discriminate G and P genotypes, respectively. VP7 and VP4 
play major roles in eliciting the production of neutralizing 
antibodies in the host immune response to rotavirus infection 
[4]. At least 27 G and 37 P genotypes have been reported to date 
[5,6]. Of these genotypes, epidemiological studies have shown that 
five G (G1-G4 and G9) and three P (P[4], P[6], and P[8]) 
genotypes are the most prevalent genotypes associated with global 
human rotavirus infections [7] . The most common G/P genotypes 
are G1P[8], G2P[4], G3P[8], G4P[8], and G9P[B] [7-9], and 
emergence of the G12 genotype in combination with P[8] or P[6] 
was recently reported [10]. A rare G/P combination, G3P[9], was 
frrst detected in humans in Japan in 1982 (AU-1 strain) [1 1,12,13]. 



Such G3P[9] strain was detected in a cat in Australia in 1984 
(Cat2 strain) and has since been commonly detected in cats [14— 
1 7] . Recently, reassortants of feline/human or feline/ canine 
G3P[9] rotaviruses have been reported in humans, which has 
serious implications for future strategies aimed at rotavirus vaccine 
development [18-21]. 

In South Korea, rotavirus is the most common viral agent 
causing acute diarrhea and gastroenteritis in young children, 
especially during the winter season. The distributions of the G and 
P genotypes show temporal and geographical fluctuations [22]. 
The major human G types are G1-G4 and G9, combined with the 
P types P[4], P[6], and P[8] [23,24]. Furthermore, unusual G and 
P combination genotypes, including G3P[3], G3P[9], G3P[10], 
G8P[8], G11P[4], G11P[25], and G12P[6], have been isolated 
from South Korea [25]. From 2003 to 2005, the G3 genotypes 
were predominant [26]. Although the P[9] genotype has been 
detected in some rotavirus surveillance studies in South Korea, no 
significant emergence of P[9] genotypes has been observed in the 
Korean population [22]; however, no studies have specifically 
focused on this genotype. Therefore, further analysis of the 
circulating Korean genotype P[9] is needed to understand its 
origin, genetic variation, and potential impact on the human 
population. 

Recently, a complete genome classification based on a 
nucleotide sequence analysis was defined according to established 
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nucleotide percent cutoff values for rotavirus [27]. According to 
this classification, the acronyms Gx-P[x]-Ix-Rx-Cx-Mx-Ax-Nx- 
Tx-Ex-Hx are used to classify VP7-VP4-VP6-VP1-VP2-VP3- 

NSPl-NSP2-NSP3-NSP4-NSP5/6-encoding gene segments, re- 
spectively [11,12]. In the current study, the G3P[9] rotavirus 
strain, RVA/Human-wt/KOR/CAU 1 2-2-5 1 /20 1 2/G3P[9] 
(CAU12-2-51), was isolated using a cell culture system. The 
complete genomic sequence of this G3P[9] virus was determined 
after two cell culture passages and analyzed in order to provide a 
better understanding of the origin and genetic composition of this 
viral strain. 

Materials and Methods 

Ethics statement 

The stool samples used in this study were collected and analyzed 
according to methods detailed in our study protocol (number 
#2011-10-06), which was approved by the Human Subjects 
Institutional Review Board (IRB) of Chung-Ang University 
College of Medicine, Seoul, Korea. Written informed consent 
was obtained from all participants in this study. For the children 
enrolled in this study, written informed consent was obtained from 
their parents. The provision of informed consent also included 
permission to use the data for future research purposes under 
conditions of anonymity. 

Patient history and rotavirus identification 

The stool specimen was collected from a 9-year-old girl who was 
hospitalized for severe gastroenteritis in August 2012 at Chung- 
Ang University Hospital, Seoul, South Korea. The patient was 
hospitalized for three days with episodes of watery diarrhea (three 
times/ day), pain in the right upper quadrant of the abdomen, and 
vomiting (nonprojectUe and nonbilious, containing consumed 
food). At that time, the patient had a hemoglobin level of 13.8 g/ 
^lL, a white blood cell (WBC) count of 8,058 WBCs/^iL, and a 
platelet count of 236,000 cells/|tL. The stool WBC count was 0- 
1 /high-power field. The patient had not received any prior 
rotavirus vaccines and had not had any contact with animals (e.g., 
cats, dogs, or cows) at least 30 days before hospitalization. Human 
group A rotavirus strain was then detected from the stool sample 
by both enzyme-linked immunosorbent assay (ELISA) and reverse 
transcription-polymerase chain reaction (RT-PCR). This strain 
was named CAU 12-2-51. The stool culture produced no growth 
for pathogenic bacteria, such as Salmonella and Shigella, which 
commonly cause diarrhea. 

Rotavirus isolation 

MA 104 cells were acquired from the Korean Cell Line Bank 
(Seoul, South Korea) and grown in minimum essential medium- 
alpha (MEM-ot; Gibco BRL, Grand Island, NY, USA) containing 
5% fetal bovine serum (FBS; Gibco BRL) and 0.1% gentamicin 
(Gentamicin Reagent Solution; Gibco BRL) at 37°C in the 
presence of 5% CO2. A G3P[9]-positive stool sample was diluted 
10-fold in phosphate-buffered saline (PBS; pH 7.4) and clarified by 
centrifugation at 10,000 for 10 min. The supernatant was 
filtered using a 0.45-pm sterile syringe filter (Corning Costar, 
Corning, NY, USA), treated with 10 |J.g/mL Trypsin 250 (Becton 
Dickinson, Sparks, MD, USA) for 30 min at 37°C. The 
supernatant was then inoculated onto MAI 04 cells in glass tubes 
and incubated with MEM-ot in the presence of trypsin (5 ng/mL), 
with constant rotation during incubation. Cells were harvested 5-7 
days after infection and subsequently passaged with MA 104 cells 
until a cytopathic effect (CPE) was achieved. The rotavirus 



remaining in the culture fluid aftc-r the final passage was examined 
by immunofluorescence and RT-PCR. 

RNA extraction and RT-PCR 

Viral double-stranded RNA was extracted using TRIzol reagent 
(Gibco BRL Life Technologies) according to the manufacturer's 
instructions. Extracted RNA was resuspended in RNase-free water 
and stored at — 80°C until use. The extracted RNA was then 
subjected to one-step RT-PCR using primer sets for the fuU 
amplification of all 11 gene segments (Table SI) [12,28,29,30,31]. 

Nucleotide sequencing and sequence analysis 

AU amplified PGR products were purified using a QIAquick 
PGR Purification Kit (Qiagen, Westburg, Germany) and then 
sequenced using a BigDye Terminator Cycle Sequencing Kit and 
an automated DNA sequencer (Model 3730; Ajjpiic-d Biosystems, 
Foster City, CA, USA). The resulting open reading frame (ORE) 
of the genes and deduced amino acid (aa) sequences were aligned 
using the CLUSTAL_X 1.81 program [28] and Lasergene 
software (DNASTAR, Madison, WI, USA). Data were compared 
with the corresponding rotavirus serjuences from the National 
Center for Biotechnology Information GenBank database. The 
nucleotide (nt) sequences obtained in this study were deposited in 
GenBank under the accession numbers KJ187594— KJ187604. 

Phylogenetic analysis 

The nt sequences of the 11 segments of the CAU 12-2-51 strain 
were compared with representative rotavirus sequences available 
from the GenBank database. Phylogenetic trees were constructed 
using neighbor-joining algorithms [25] from the PHYLIP suite [29] 
and the Kimura two-parameter model using MEGA5.03 software 
[30,31]. Evolutionary distances for the neighbor-joining analysis 
were based on the model described by Jukes et al. [32]. Tree 
topology was evaluated using the bootstrap resampling method with 
1 000 replicates of the neighbor-joining dataset with the SEQBOOT 
and CONSENSE programs from the PHYLIP suite. 

Results 

Rotavirus genotyping and isolation 

The VP7 and VP4 nt sequences of the CAU 12-2-51 strain 
obtained from the first-round PGR products were identical to the 
sequences of G3 and P[9] genotypes available in the GenBank 
database. CPEs were observed after two culture passages when 
MA 104 cells exhibited obscure borders, fusion, rounding, and 
detachment from the wails of the tubes (data not shown). 

Gene sequence comparisons and phylogenetic analysis 

The complete full-length nt and deduced aa sequences of the 
CAU 12-2-51 strain were obtained. For reference, only the strains 
that contained [:omplete or at least 50'/^) sequence availability of 
ORFs, as suggested by Matthijnssens et al. [27], were used. 

FuU genome-based classification. The CAU 12-2-51 strain 
was found to possess the G3-P[9]-I2-R2-C2-M2-A3-N2-T3-E3- 
H3 genotype constellation for die VP7-VP4-VP6-VP1-VP2-VP3- 
VP4-NSP1-NSP2-NSP3-NSP4-NSP5/6 genes, respectivefy 
(Table 1, Figure 1). The VP7, VP4, .NSPl, .mP3, j\SP4, and 
NSP5 genes were related to the human/feline AU-l-like genotype 
and Australian feline strains. The remaining five genes, namely, 
VP6, VP1-VP3, and JVSP2, were related to die human/bovine DS- 
1-like and bovine-like G6 human strains. 

VP4 gene. Phylogenetic analysis of the P[9] genotypes 
revealed five distinct clusters, in which the VP4/P[9] genotype of 
the prototype strain AU-1 was assigned to cluster a. Phylogenetic 
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Figure 1 . Phylogeny of the 1 1 genomic segments of tlie CAU1 2-2-51 strain. The different genotypes within each segment are indicated. 
Numbers at the nodes indicate the level of bootstrap support (%) based on neighbor-joining analysis of 1 000 resampled datasets. Only values above 
70% are provided. The CAU1 2-2-51 strain is marked in bold type and indicated by filled circles. Color coded boxes are used to differentiate strains of 
origin: human (light blue), bovine (green), feline (red), human/bovine-like (yellow), human/feline-like (purple), feline/bovine-like (blue), canine/ 
bovine-like (orange), guanaco/feline-like (dark blue), and antelope/bovine-like (light green). Genotype origins are designated as follows: G, 
glycosylated; P, protease sensitive; I, intermediate capsid shell; R, RNA-dependent RNA polymerase; C, core shell protein; M, methyltransferase; A, 
interferon antagonist; N, NTPase; T, translation enhancer; E, enterotoxin; H, phosphoprotein. 
doi:l 0.1 371/journal.pone.0097127.g001 



analysis showed that the VP4/P[9] genotype of the CAU12-2-51 
strain was classified into cluster c, together with the recently 
identified human and feline P[9] genotypes, including the feline 
P[9] strain isolated in Italy (BA222) and human/feline-like strains 
isolated in France (RI320), Japan (KF17), and Russia (0211) 
(Figure la). Interestingly, VP4/P[9] of the CAU12-2-51 strain 
showed the highest nt and aa sequence identities to the G3P[9] 
feline/human reassortant BA222 strain (97.9% and 99.2%, 
respectively). On the other hand, VP4/P[9] of the CAU 12-2-51 
strain showed lower sequence identity to VP4/P[9] of the prototype 
strain AU-1 (95.4% at the nt level and 97.3% at the aa level). 

VP7 gene. Phylogenetic analysis revealed that the VP7/G3 
genotype of the CAU 12-2-51 strain was clustered into lineage 3, 
together with recent human G3P[9] strains isolated from Thailand 
(CU365-KK, CMH120, and CMH134) and China (L621), 
showing the highest nt and aa sequence identities (>96% and > 
97.5%, respectively; Figure lb). In contrast, the VP7/G3 genotype 
of the CAU 12-2-51 strain showed lower nt and aa sequence 
identities (90.2% and 94.1%, respectively) to the G3 prototype 
strain AU-1 (Tabic 2). 

VPl, VP2, VPS, VP6, and NSP2. Matthijnssens et al. 
demonstrated that human rotavirus strains bcdong to the DS-1- 
like genogroup (G2-P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2), shar- 
ing the VP6, VPl, VP2, VP3, NSP2, and NSP4 genotypes (I2-R2- 
C2-M2-N2-E2) with those of bovine rotavirus strains [12]. In this 
study, genotypes of VP6-VP1-VP2-VP3-NSP2 genes of the 
CAU 12-2-51 strain were classified within the DS-l-like genogroup 
(I2-R2-C2-M2-N2), commonly seen in human/bovine-Kke and 
bovine G6 rotavirus strains (Figure Ic-f and Ih). However, the 
sequence identities of these genes to DS-1 were relatively low, i.e., 
just above the cutoff values, showing nt (aa) sequence identities of 
86.2% (98.4%), 8(i.3% (97.4%), 84.8% (96.9%), 83.6% (90.9%), 
and 86.2% (93.6%) for VP6, VPl, VP2, VPS, and MP2 genes, 
respectively (Table 2). The VP6, VPl, VP2, VPS, and NSP2 genes 
of the CAU 12-2-51 strain showed the highest nt (aa) sequence 
identities to sequences of human/bovine-like rotavirus strains 
isolated in Hungary (strain Hun5 for VPl and VP2), Japan (strain 
KF17 for VPS), Russia (strain RUS-Nov-K2 for VP6), and Russia 
(strain 0211 ior: NSP2) at 98.5% (99.7%), 97.4% (99.1%), 97.2% 
(99.4%), 97.8% (97.2%), and 99.4% (99.6%), respectively 
(Table 2). The VP6 VPl gene segments of the CAU 12-2-51 
strain clustered closely with strains originating from human/ 
bovine-like rotavirus strains (BP1879, Hun5, KF17, PAH136, 
PAIll-96, NovlO-N507, and RUS-Nov06-K2), a feline /bovine- 
Uke G3P[9] strain (BA222), and a human/bovine-porcine like 
G12P[4] strain (L26; Figure Ic, Id, and Table S2) [19,33-37]. In 
contrast, VP2, VPS, and NSP2 gene segments were found to cluster 
together with most of the known rotavirus strains originating from 
human/bovine-like G6P[14] rotaviruses (B1711, Hun5, PAH136, 
PAI58, and KF17), a feUne/bovine-like G3P[9] rotavirus (BA222), 
and bovine rotaviruses (NCDV, RF, WC3, 1063, and 1065; 
Figure le. If, Ih, and Table S2) [12,19,37,39,40,42,43,44,45,46]. 

NSPl, NSP3, NSP4, and NSP5. Based on RNA-RNA 
hybridization, the human AU-l-like genogroup (G3-P[9]-I3-R3- 
C3-M3-A3-N3-T3-E3-H3) is beheved to have a close evolutionary 
relationship with canine and feline rotavirus strains [12]. Our 



study indicated that the NSPl, MSPS, .NSP4, and .NSP5 genes of the 
CAU 12-2-51 strain were classified with the genotypes of strains 
with AU-l-like genotypes, showing nt (aa) sequence identities of 
93.1% (93.1%), 96.9% (97.7%), 96.2% (97.1%), and 93.1% 
(92.9%), respectively (Figure Ig and li-k). On the other hand, the 
MPl, NSP4, and NSP5 genes of the CAU12-2-51 strain showed 
the highest sequence identities with the G6P[9] genotype of the 
human/feline-like KF17 strain, with nt (aa) identities of 97.2% 
(96.7%), 97.7% (97.1%), and 98.8% (97.9%), respectively. 
Additionally, the NSP3 gene showed the highest nt (aa) sequence 
identity to tiie AU-1 strain at 96.9% (97.7%) (Table 2). 
Phylogenetic analysis indicated tiiat the NSP1-NSP3-NSP4- 
NSP5 gene segments found in A3-T3-E3-H3 genotypes clustered 
closely with the strains originating from human/feline-like 
rotaviruses (K8, KF17, PAH136, PAI58, and RV10109), a 
guanaco/feline-like rotavirus (Chubut), and feline rotaviruses 
(Cat2 and BA222; Figure Ig, li-lk, and Table S2) 
[18,19,33,35,36,40,41]. 

Discussion 

Rotavirus infection causes acute diarrhea and is the most 
common cause of gastroenteritis affecting infants and young 
children in South Korea. The human rotavirus G1-G4 genotypes 
account for almost all rotavirus infections circulating in South 
Korea in recent years. Gl was the predominant genotype before 
2000 and between 2004 and 2009; G2 and G3 were more 
prevalent during the seasons of 2000-2001 and 2003-2005, 
respectively; and G4 was predominant during the seasons of 2002- 
2003 and 2007-2008 [42]. G3 rotavirus genotypes have a broad 
host range and they have been discovered in several host species. 
G3 has been shown to be the dominant genotype and has spread 
throughout many Asian countries, including Japan [43] , China 
[44], Vietnam [45], and Hong Kong [46]. Recently, the 
incidences of the G2, G3, and G4 genotypes decreased, while 
the incidence of the Gl genotype increased, suggesting that the 
predominant genotypes of rotaviruses in South Korea can change 
rapidly within a short period of time [33] . 

The P[9] genotype is frequentiy detected in cats [47], and the 
first P[9] human rotavirus, the K8 strain with Gl specificity, was 
isolated from a 14-year-old boy in Hokkaido, Japan in 1977 [48]. 
Although P[9] genotypes have been detected sporadically in most 
of the rotavirus surveillance studies conducted in South Korea, 
there has been no significant emergence of P[9] genotypes in the 
Korean population. The relative global frequency of the P[9] 
genotypes represents less than 2.5% of total rotavirus infections, 
and the majority of these genotypes carry Gl, G2, and G3 
specificity, except for a few isolates that have been found to have 
G4, G6, G9, and G12 specificity [20]. The present analysis of 
VP4/P[9] genes revealed a genetic relationship between the 
CAU 12-2-51 strain and the recently identified P[9] strains, 
demonstrating that this strain had the highest identity to the 
feline rotavirus strain BA222. 

Recent reports established a rotavirus classification system in 
which each gene of the rotavirus is assigned a particular genotype 
based on designated nucleotide identity cutoff values. Full genome 
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analysis also provided an excellent platform for determining 
interspecies evolutionary relationships. To date, complete genome 
characterizations have assigned the G3P[9] strains to the genotype 
constellations G3-P[9]-I3-R3-C3-M3-A3-N3-T3-E3-H3 [12], G3- 
P[9]-I2-R2-C2-M2-A3-N1-T6-E2-H3 [36], G3-P[9]-I2-R2-C2- 
M2-A3-N2-T6-E2-H3 [36], G3-P[9]-I3-R3-C3-M3-A3-N3-T3- 
E3-H6 [7,21], G3-P[9]-I2-R2-C2-M2-A3-N1-T3-E2-H3 [9], 
G3-P[9]-I3-R3-C2-M3-A3-N1-T6-E3-H3 [18], and G3-P[9]-I2- 
R2-C2-M2-A3-N2-T1-E2-H3 [49]. The genotype of the human 
CAU12-2-51 strain exhibited a new G3-P[9]-I2-R2-C2-M2-A3- 
N2-T3-H3 genotype constellation. 

Elucidation of the complete genomic sequence of the C AU 12-2- 
51 strain revealed that the genes encoding VP7, NSPl, NSP3, 
NSP4, and NSP5 were related to the human AU-l-like genotype, 
together with the related human/feline-like and feline rotavirus 
strains. The remaining five genes encoding VP6, VP1-VP3, and 
NSP2 were related to the human DS-l-like genotype and the 
related human/bovine-like and bovine rotavirus strains. These 
results were in agreement with the hypothesis that genotypes of 
human rotavirus strains belonging to the AU-l-like and DS-l-like 
genogroups had a close evolutionary relationship with feline/ 
canine and bovine rotavirus strains, respectively [12]. However, it 
is currently difficult to determine the specific origins of some of 
these genotypes because of the unavailability of sequence data for 
animal rotavirus strains, including feline and bovine rotaviruses. 
Therefore, additional sequence data for these strains are needed to 
understand the true origins of these viruses [12]. 

There are two approved vaccines for rotavirus, which were 
developed based on diflFerent strategies. Rotarrx (GlaxoSmithKline 
Biologii:als, RLxcnsart, Belgium) is a single, live-attenuated human 
G1P[8] rotavirus strain [50], while RotaTeq (Merck and Co. Inc., 
West Point, PA) is a pentavalent, live bovine-human reassortant 
vaccine containing the Gl, G2, G3, G4, and P[8] genotypes [51]. 
RotaTeq and Rotarrx were introduced in South Korea in 
September 2007 and July 2008, respectively. The patient in this 
study had not received any rotavirus vaccines and showed severe 
symptom of diarrhea. Comparison of the aa sequence of the VP4 
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